Fatigue in disordered media 
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We obtain the Paris law of fatigue crack propagation in a disordered solid using a fuse network model where 
the accumulated damage in each resistor increases with time as a power law of the local current amplitude. 
When a resistor reaches its fatigue threshold, it burns irreversibly. Over time, this drives cracks to grow until 
the system is fractured in two parts. We study the relation between the macroscopic exponent of the crack 
growth rate - entering the phenomenological Paris law - and the microscopic damage-accumulation exponent, 
y, under the influence of disorder. The way the jumps of the growing crack, Aa, and the waiting-time between 
successive breaks, At, depend on the type of material, via y, are also investigated. We find that the averages of 
these quantities, (Aa) and (At) / (t r ), scale as power laws of the crack length a, (Aa) °= a a and (At) / (t r ) oc a~$, 
where (t r ) is the average rupture time. Strikingly, our results show, for small values of y, a decrease in the 
exponent of the Paris law in comparison with the homogeneous case, leading to an increase in the lifetime of 
breaking materials. For the particular case of y = 0, when fatigue is exclusively ruled by disorder, an analytical 
treatment confirms the results obtained by simulation. 



Fatigue failure is relevant in most everyday-life situations. 
For instance, asphalt pavement has its lifetime decreased by 
cyclic traffic loading, and aircraft fuselage has its lifespan de- 
termined by the number of pressurization cycles [!1]. Fatigue 
is the largest cause of engineering failure, comprising an esti- 
mated 90% of all metallic failures J2l- Furthermore, due to its 
non-linear and critical behavior, it is mostly catastrophic and 
insidious, occurring very suddenly and without warning if the 
crack is not noticed, since it is brittle-like even in normally 
ductile materials J^IQ]. So far, theoretical understanding of 
crack-growth phenomena is far from complete, and fatigue- 
life prediction remains essentially empirical [5]. 

Generally speaking, failure in stressed materials is the re- 
sult of accumulated microscopic damage due to cyclic or con- 
stant external loads. For cyclic (constant) load one speaks of 
fatigue (creep) failure Ja,|7|]. In a mesoscopic description, 
when the accumulated microscopic damage reaches a local 
threshold, cracks are initiated and grow through the material, 
leading to macroscopic fracture. Under fatigue, local dam- 
age increases and cracks grow even for external loads below 
the limit strength of the material. Crack growth in this sub- 
critical regime, for an intermediate crack length a, is usually 
described by the empirical Paris law |@], 
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which relates the crack velocity da / dt with the stress-intensity 
amplitude AK Aa^/a, for a sample subject to an external 
stress amplitude Aa. The constants C and m depend on the 
material and on the loading conditions QJJ] . 

Due to the intrinsic characteristics of random microscopic 
defects in fracturing materials |9 f4l6ll . statistical models such 
as the random spring model 11711 . the random beam model 
dHUH], an d the random fuse model Ji, 1^-21] have been suc- 
cessfully applied to study fracture. The role of disorder in 
these systems is to induce a transition from brittle to duc- 
tile behavior Ji, 22-25]. As already mentioned, the presence 



of fatigue can also lead to brittleness, although the physical 
mechanism responsible for this behavior so far remains un- 
known [2]. Vieira et al. [26;] have shown analytically and 
confirmed numerically the relation between the Paris-law ex- 
ponent m and the microscopic damage-accumulation param- 
eter y in homogeneous materials. In the same study, first at- 
tempts have been made to analyze the stability of the obtained 
results towards the introduction of disorder. The aim of the 
present work is to provide the full picture, using fuse-model 
simulations and analytical arguments, of the role of disorder in 
fatigue crack-growth phenomena, including the Paris as well 
as the catastrophic (post-Paris) regime. 

In our implementation, the random fuse model consists of 
resistors (with the same conductance) located on the bonds of 
a tilted square lattice, mimicking springs in the more compli- 
cated case of fracture mechanics 1141 12311 . In this scalar version 
of elasticity, the current (voltage) in the electrical system is 
analogous to the magnitude of the stress (strain) in the me- 
chanical system. In order to simulate heterogeneous micro- 
scopic defects commonly found in real materials, we assign 
to each resistor a random fatigue threshold F th , chosen from 
a uniform probability distribution between 1 — b and 1 + b, 
where b gauges the strength of disorder. Initially, the sys- 
tem has a notch, formed by removing two contiguous bonds, 
placed symmetrically at the center of the system, with peri- 
odic boundary conditions along the notch direction. A global 
transverse current is imposed across the system from top to 
bottom and the local currents are numerically calculated by 
applying Kirchhoff 's law to each node 12711 (nodes being the 
points where four resistors meet). We then compute, for each 
resistor i, the fatigue damage history up to time t by using a 
power-law damage-accumulation function, 
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where /, (?) is the current in resistor i at time t, and the stress- 
amplification exponent y is a phenomenological microscopic 
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FIG. 1 . (Color online) Snapshots of the current distribution over the 
lattice for y = 3 (above) and y = (below) with b = 0.5 and L = 256. 
Initially, the system has a notch at its center, which grows horizon- 
tally, fracturing the medium in two parts. Currents increase as colors 
change from yellow, over violet, to green. The main crack, the one 
connected to the notch, of length a, is colored black, while secondary 
cracks are shown in red. 



parameter, dependent on the material. Below the fatigue 
threshold, a resistor follows Ohm's law, but once fatigue dam- 
age reaches the corresponding threshold, /,■ = F\ h , the resistor 
burns irreversibly, like an electrical fuse. The system here is 
operated in the quasi-static limit, i.e., the period of the current 
cycle is much larger than the time required for the rearrange- 
ment of the current distribution following the burning of a 
fuse, but much smaller than the time between two consecutive 
burning events. This is equivalent to assuming a well-defined 
current amplitude A/, in each resistor during that time. There- 
fore, the damage increment at resistor i between two consec- 
utive burning events occurring at times t and t + At is simply 
given by 



4/i(Af) = AAt [Ali 



(3) 



where A/,- (f ) is the current amplitude in resistor i between 
events and the constant A sets the time scale so that the dam- 
age fi becomes dimensionless. 

Only the L resistors in the midline of the system, which in- 
cludes the initial notch, are allowed to burn (see Fig.[T|). Thus, 
the results reported here are valid as long as the crack tends 
to propagate perpendicularly to the direction of the external 
load. After a resistor burns, the current distribution is rear- 
ranged, causing an increase in the current around crack ends. 
We allow only one resistor to burn at a time and this burn- 
ing process continues until a macroscopic fracture crosses the 
system. Due to the presence of disorder, the main crack, start- 
ing from the notch, can grow either by one resistor at a time 
or by coalescence with secondary cracks initiated ahead of the 



notch. In the latter case, it advances several resistors instanta- 
neously, stopping at the next resistor which has not reached 
its fatigue threshold. This procedure of following a single 
one-dimensional crack has been used experimentally to study 
growing cracks in thermally activated stressed systems 12811 . 
As Fig.[T]shows, increasing the value of y greatly reduces the 
occurrence of secondary cracks. 

In the special case of y = 0, the damage-accumulation rate 
is equal for all resistors, and independent on the current am- 
plitude, and the breaking sequence is determined exclusively 
by the random distribution of fatigue thresholds. This partic- 
ular limit of y is especially important as it enhances the role 
of disorder in fatigue failure. The process is then similar to 
one-dimensional percolation, and analytical equations can be 
developed to confirm our results from numerical simulations. 

For y = 0, one needs only to distinguish between the ho- 
mogeneous (b = 0) and finite disorder (b ^ 0) cases, since 
systems with different values of b > can be mapped on each 
other by rescaling time. Thus, we consider a one-dimensional 
chain of resistors whose fatigue thresholds are chosen accord- 
ing to a uniform probability distribution, P (F th ), in the inter- 
val [0,1]. We define the time unit by imposing that the accu- 
mulated damage at each resistor at time t is /(f) = t, which 
is equivalent to measuring time in units of the rupture time. 
Thus, the probability that a given resistor is already burnt at 
time t is simply t. Given that the main crack has length a, 
having started from a pair of adjacent burnt resistors at time 
t = 0, the probability that it advances Aa resistors between 
times t and t + dt, after having waited a time between At and 
At + d(At) since it last advanced, is P(t,At,Aa\a) dtd(At), 
with 



P{t,At,Aa\a) = (a - 2) (a - I) a (t - Atf~ 3 (1 - 1 )f A ""' 

x [(1 -t) (1 - 8 A „, L _„_i) + aAfl^a-i] • (4) 

In order to arrive at the above conditional probability, let 
us assume that the crack advances to the right, and notice 
that (t - At)"' 3 (1 - f)d(Af) is the probability that a - 3 re- 
sistors burned before time t — At, one resistor burned between 
t — At and t — At + d(At) (while two resistors form the ini- 
tial notch), and the resistor at the left end of the crack is still 
intact at time t. The factor t ~ l dt corresponds to the proba- 
bility that the crack advances Aa resistors between times t and 
t + dt, while the terms in square brackets differentiate a crack 
that advances less than L — a — 1 resistors, stopping at yet an- 
other intact resistor, from a crack that goes around the system 
(due to periodic boundary conditions). Of course the process 
is symmetric for a crack advancing to the left, and the fac- 
tor (a — 2) (a — 1) a ensures normalization of the conditional 
probability. 

The marginal probabilities P(Aa\a) andP(Af |a), which cor- 
respond to the distributions of crack jumps and waiting times 
between successive jumps, for a given crack length a, can be 
readily obtained from P(t,At,Aa\a) by integrating over the 




a/L 



FIG. 2. (Color online) Main plot: Rescaled average jump of the 
main crack as a function of its length a, for various system sizes, two 
values of the damage accumulation exponent, y = and y = 3, and 
disorder strength b = 0.5. The results for y = match the analytical 
prediction (solid line) from Eq. Q. The jump of a crack depends 
on its length as a power-law, (Aa) fl oc a a . The error bars are smaller 
than the symbols. The inset shows the exponent a versus y for dif- 
ferent degrees of disorder. In the ordered case, b — 0, the exponent a 
vanishes for any value of y since the crack advances without coales- 
cence. 



~ i — i i i i 1 1 1 

y = 3 



o 


L = 512 


A 


L = 256 




L = 128 


* 


L = 64 




Eq. (8) 



ex 




a/L 



FIG. 3. (Color online) Averaged waiting time between successive 
breaks of the main crack, normalized by the averaged rupture time, 
(t r ) , as a function of the length a, collapsed for different system sizes, 
for y = and y = 3, with b = 0.5. The results of y = match the 
analytical prediction (solid line) from Eq. {8}. The waiting-time of a 
crack depends on its length as a power-law, (Ar) fl / (t r ) °= aT$ . The 
error bars are smaller than the symbols. The inset shows the exponent 
(3 versus y for different degrees of disorder. 



remaining variables, and are given by 

p / * I \ 2(a-l)a(l-8 AuX - fl -i) 

P (Aa \a) = — — — 

(a + Afl-2)(a + Afl- l)(a + Aa) 

(g-l)g s rz\ 

+ (L-3)(L-2) W -"- 1 (5) 

and 

P(At\a) =a{\ -At)"' 1 . (6) 

From these, we can calculate, as functions of the crack length 
a, the average crack jump and the average waiting time be- 
tween consecutive jumps, 

(Aa) a = £ AaP(Aa\a)=a(l-^-) (7) 

Aa=l V L L J 

and 

(At) a = f l d(At)AtP(At\a) = -^—. (8) 
Jo a + 1 

Notice that (Aa) a grows linearly with a for 1 C fl C t, but 
decreases rapidly as a approaches L, since Aa cannot be larger 
than L — a — 1. On the other hand, (At) a decreases as l/a for 
o> 1. 

In the presence of disorder, there are various possible defi- 
nitions for the crack velocity. For instance, we can determine 
the average time dependence of the crack length, (a(f)), and 
calculate its derivative. Alternatively, for a given crack length, 
we can estimate crack velocity in terms of averages as 

/Aa\ (Aa). I AA" 1 

W fl ' (A/) fl ' \Aa/ a • 



However, and not surprisingly, it is possible to show that, for 
y = 0, all these definitions yield crack velocities which scale 
with the crack length a, for 1 <C a L, as 

v~a\ (9) 

leading, as defined in Eq. (Q~|), to a Paris exponent m = 4. In 
our simulations, we estimate the crack velocity, da/dt, by a 
procedure designed to minimize statistical fluctuations. Pre- 
cisely, for a given disorder realization, we estimate the time at 
which the crack had a given length. (For disorder realizations 
in which the crack happened to be temporarily trapped at that 
particular length, we chose the time at which the crack ad- 
vanced from that length.) Then, we calculate the average time 
at each crack length over all disorder realizations. We take av- 
erages over up to 50 000 realizations, depending on the set of 
parameters, y, b, and L. Finally, we numerically differentiate 
the corresponding curve with respect to time, to estimate the 
crack velocity. 

We now discuss the results obtained by simulating the fuse- 
network model with y > 0, and compare the case of y = to 
the above analytical expressions. The jumps and waiting times 
of the main crack, as functions of the length a, are analyzed 
via finite-size scaling in Figs. [2] and [3] respectively. The re- 
sults show that both (Aa) a and (At ) a / (t r ) scale as power laws 
of the crack length a, (Aa) a °< a a and (At ) a / (t r ) °= a~$, where 
(At) a is normalized by the average rupture time (t r ). The be- 
haviors of the exponents a and p as functions of the stress- 
amplification exponent y, for different disorder strengths b, 
are shown in the insets of Figs. |2] and [5J respectively. For 
y = and b>0, both a and p are unity, in agreement with the 
analytical results. For the ordered case (b — 0), a vanishes for 
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FIG. 4. (Color online) Rescaled crack growth rate of the main crack 
as a function of rescaled length a/L, showing the Paris exponent, m, 
for y = and y = 3, with b = 0.5. For better visualization, curves for 
y= 3 are shifted to the right by 20%. 




FIG. 5. (Color online) The Paris exponent m versus the stress- 
amplification exponent y for different disorder strengths b. The 
solid lines represent the results found for the homogeneous case, 
m = 6 — 2y, for y < 2, and m = y, for y > 2 |2(J. For the one- 
dimensional case with y = 0, m = 4 for any finite value of b. 



all values of y, since the main crack always advances two re- 
sistors at a time, and no secondary cracks appear. In this limit, 
P is one half of the Paris exponent m. 

The growth rate of the main crack is depicted in Fig. |4j 
which shows the rescaled crack velocity for different system 
sizes, disorder strength b = 0.5, and y = and 3. The values 
of the Paris exponent are chosen so as to yield the best data 
collapse of the curves corresponding to different system sizes 
for a given value of y. Usual materials present, after the Paris 
regime, a faster growth rate preceding the final rupture. How- 
ever, due to boundary effects, we observe, for any finite disor- 
der strength (b > 0) and for y close to zero, a slight decrease in 
the growth rate as the system approaches catastrophic failure. 

The dependence of the macroscopic Paris exponent m on 
the microscopic fatigue exponent y, for different degrees of 
disorder, is shown in Fig. [5] together with the results found 
for the homogeneous case, m = 6 — 2y, for y < 2, and m = y, 
for y > 2 12611 . Numerical estimates extracted from fuse-model 
simulations for the homogeneous case are also shown, high- 
lighting the existence of logarithmic corrections to finite-size 
scaling, which plague the results for values of y around the 
critical value y c = 2. In the presence of disorder (b > 0), m 
converges to 4 as y approaches zero, as predicted by the an- 
alytical results, while the behavior is similar to the ordered 
case for y > 2, except for very strong disorder (b ~ 1). Inter- 
estingly, any finite disorder leads to a decrease in the Paris ex- 
ponent, and thus to a potential increase in the fatigue lifetime, 
for y < 1 . On the other hand, in agreement with a stability 
analysis of the effects of disorder on the homogeneous system 
02611 . the system is essentially insensitive to small disorder for 
y>2. 

In conclusion, we have investigated both analytically 
and by fuse-model simulations the behavior of fatigue fail- 
ure in heterogeneous materials, characterized by a damage- 



accumulation law governed by a stress amplification exponent 
y. We have determined how the empirical, macroscopic Paris 
exponent depends on the microscopic parameters of fatigue 
and disorder, represented by random local fatigue thresholds. 
For small values of y, corresponding to materials in which 
damage accumulation depends only weakly on the local stress 
amplitude, the Paris exponent approaches m = 4 (obtained ex- 
actly for y = 0) if any disorder is present. On the other hand, 
fatigue failure in materials for which y > 2 continues to be 
characterized by a Paris exponent m — y as long as disorder 
is not sufficiently strong. A considerable challenge for fu- 
ture research lies in obtaining the relation between atomistic 
processes and the coarse-grained description of damage accu- 
mulation by means of a stress-amplification exponent y. 
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